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ABSTRACT 

Periodicity search in gamma-ray data is usually difficult because of the small number 
of detected photons. A periodicity in the timing signal at other energy bands from 
the counterpart to the gamma-ray source may help to establish the periodicity in the 
gamma-ray emission and strengthen the identification of the source in different energy 
bands. It may, however, still be difficult to find the period directly from X-ray data 
because of limited exposure. We developed a procedure, by cross-checking two X-ray 
data sets, to find candidate periods for X-ray sources which are possible counterparts to 
gamma-ray pulsar candidates. Here we report the results of this method obtained with 
all the currently available X-ray data of 8 X-ray sources. Some tempting periodicity 
features were found. Those candidate periods can serve as the target periods for future 
search when new data become available so that a blind search with a huge number of 
trials can be avoided. 

Key words: pulsars: general - stars: neutron - gamma-rays : observations - X-rays: 
general 



1 INTRODUCTION 

The electrodynamics in the magnetosphere of a neutron star 
has not yet been well understood, although fascinating phe- 
nomena as revealed in emissions across the whole electro- 
magetic spectrum from neutron stars have been observed 
for forty years. Among the 2000 or so radio pulsars cata- 
loged to date, only 7 are detected at photon energies higher 
than 100 M eV and 3 more a re detected at a lower confi- 
dence level (|Thompsonll2004l) . These 7 include Geminga, 
which do es not emit radio waves or emits only at a very 
low level l|Kuz'min fc Losovskiilll999l ). The effort to under- 
stand gamma-ray emissions from puls ars has resulted in two 
types of models: the slot-gap models dMuslimov fc Hardind 
20041) and the outer-gap models (e.g. [ Cheng et al.l (|2000l ) ; 
Hirotanil (l2007l ); lTakata fe Chanel l|2007l )). These models, al- 
though different in the assumed emission site and mecha- 
nisms, have been successful in many aspects in explaining 
the pulse profiles and spectra. However, they do give dif- 
ferent estimates of the fraction of radio-quiet gamma-ray 
pulsars, like Geming a, among all the gamma-ray pulsars 
(|Harding et al.ll2007l ). A gamma-ray pulsar may be radio- 
quiet if its gamma-ray beam is wider than its radio beam 
or these beams are in different directions. It is also possible 
that it is intrinsically weak in radio emission because of un- 
favored spin periods and magnetic field strength. Although 
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a consistent model for radio emission from pulsars is still 
lacking, several empirica l models are available. Given a ra- 
dio emission model (e.g. IJohnston fe Weisberl (|200rj) ), dif- 
ferent models for gamma-ray pulsars may be discriminated 
by comparing their estimates of the radio-quiet gamma-ray 
pulsar fraction with observations, in addition to comparing 
their predictio n of high-energy em ission polarization char- 
acteristics (e.g. iTakata et al.l |2007)). 

The number of currently known gamma-ray pulsars is 
only 7, which is obviously too small to give a meaningful 
fraction of the radio-quiet ones. With the launch of GLAST 
and AGILE, the number of gamma-ray pulsars is expected 
to increase significantly. However, radio-quiet gamma-ray 
pulsars may still be difficult to find, mainly because of 
the limited amount of gamma-ray photons. Just like the 
case of Geminga, only after its period was discovered in 
the ROSAT data (|Halpern fe Holtl Il992l ). its pulsation in 
gamma-rays was then found. Very often it is also difficult to 
find the period in X-ray data, particularly when the data 
is sparse and a blind search over a large range of peri- 
ods is performed. Efforts to find new radio-quiet gamma- 
ray pulsars from unidentified EGRET sources in a multi- 
wave length approa c h hav e not yet yielde d defin ite results 
(e.g. iReimer et ail (|200ll ): lHalpern et ail <|2002l )). Among 
the unidentified EGRET sources, many are at low galac- 
tic latitude and may be gamma-ray pulsars. Continuing this 
effort, we conducted a periodicity search for possible X-ray 
counterparts to some bright unidentified EGRET sources, 
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using all available archival data with the method of cross- 
checking two data sets to look for candidate periods. The 
purpose of this cross-checking procedure is to find relatively 
significant signatures of possible pulsations occurring in two 
independent data sets. Candidate periods and the corre- 
sponding period time derivatives thus found are reported 
and can serve as the target periods to verify in future X-ray 
or gamma-ray data. 



2 SOURCE SELECTION AND DATA 
REDUCTON 

As an effort to understand t he nat ure of unidentified 
EGRET sources, iRoberts et alJ l)2001ah presented a catalog 
of 2-10 keV ASCA/GIS images of fields containing 30 bright 
GeV sources. These GeV sources were selected from the 
bright est ones in the catalog compiled bv lLamb fc Macomb! 
l|l997f ) using only photons with energies above 1 GeV from 
the EGRET public archives that incorporate the first 4.5 yr 
of CGRO/EGRET observations. Among these GeV sources, 
some are identified pulsars and quasars, some are possible 
counterparts to previously known objects, and the others are 
not yet identified to sources in other energy bands. Possible 
X-ray coun terparts of the s e unide ntified GeV sources were 
proposed in lRoberts et al.l l|2001al ). To search for new radio- 
quiet gamma-ray pulsars, we look for periodicity of possible 
X-ray counterparts to those unidentified GeV sources for 
which the possibility of being gamma-ray pulsars has not yet 
been ruled out. Furthermore, as explained in the next sec- 
tion, two data sets adequate for timing analysis are needed 
for the cross-checking method to find candi date periods. 6 
GeV so urces among the 30 in the catalog of IRoberts et al.l 
l|2001al ) were thus selected. One of the 6 sources (GeV J1417- 
6100) has four potential X-ray counterparts. Another one 
of the 6 sources (GeV J0008+7304) has been investigated 
with the same procedure described in this paper. The peri- 
odicity search result of its pos sible X-ray counterp art, RX 
J0007.0+7302, was reported bv lLin fc Chanel (|2005l ) and will 
not be repeated here. The 8 X-ray sources and their data 
used in this study are listed in Table [1] De tailed descrip- 
tions o f each GeV sources can be found in IRoberts et al.l 
i|2001al ). 

GeV J1417 -6100 is in the direc tion of the 'Kookaburra' 
radio complex (Ro berts et al.|[l999f) . A 68-ms radio pulsar, 
PSR J1420-6048, was found at a positi on consistent with 
AX J 1420. 1-6049 (|DAmico et all l200ll ). which is in the 
northeastern wing of the Kookaburra. The pulsation at the 
radio period in X-rays was only marginally detected with 
ASCA data (jRoberts et al.ll2001bl ). Chandra images of AX 
J1420. 1-6049 clearly show a compact source at the position 
of PSR J1420-6048, but its pulsation in X-rays at the radio 
period cannot be confirmed w ith the current Chandra data 
l|Ng. Roberts fc Romanll2005h . The identification between 
PSR J1420-6048 and AX J1420. 1-6049 is not yet conclu- 
sive. We therefore include AX J 1420. 1-6049 in our period- 
icity search. AX J1418. 7-6058 is located in the southwest- 
ern wing of the Kookaburra (the 'Rabbit'). Two compact 
sources denoted as 'Rl' and 'R2' in the Rabbit were identi- 
fied i n Chandra and XMM images (|Ng. Roberts fc Romanil 
2005). Although 'Rl' might be a background source, we still 
used it as the target in the XMM data for cross-checking 



periodicity signatures with earlier ASCA and BeppoSAX 
data, since 'Rl' is much brighter than 'R2' and in those ear- 
lier data photons from 'Rl' dominate. More recently, two 
extended TeV sources (HESS J1420-607 and HESS J1418- 
609), spatially coinciden t with the two wings o f the Kook- 
aburra, were discovered l|Aharonian et al . 2006). It is likely 
that 3EG J1420-6038 and GeV J1417-6100 are two separate 
sources associate d with HESS J1420-60 7 and HESS J1418- 
609 respectively (lAharonian et al] [2006). Another two X-ray 
sources, AX J1418. 2-6047 and AX J1418.6 -6045, denoted as 
'Src2' and 'Src3' in lRoberts et all l|2001al ). are both located 
in the 95% error contour of GeV J1417-6100 and in the 99% 
contour of 3EG J1420-6038 l|Aharonian et al.ll2006l ). 

AX J1809. 8-2332 is the brightest X-ray source within 
the 99% error box of GeV J1809-2327 (3EG J1809-2328). 
Chandra images of this region reveal a point X-ray source 
conn ected to a nonthermal X-ray/radio nebula (|Braie et al.l 
2002). XMM images show a larger, broader X-ray emission 
region with a ridge of emission along t he symmetric axis 
of the radio nebula l|Roberts et alJ [20061 ). The X-ray point 
source is very likely a pulsar moving at the tip of its asso- 
ciated pulsar wind nebula. It is therefore very much desired 
to find its pulsation period. 

GeV J1835+5921 and its possible X-ray counter- 
part AX J1836. 2+5928 have long been considered as the 
next Geminga to find. Strong evidence from Chandra 
and HST observations indicates AX J1836. 2+5928 (RX 



J1836.2+5925) is a neutron star (|Halpern et al.ll2002h. Us- 
ing th e Chandra data listed in Table 111 lHalpern et all 
l|2007l ) reported unsuccessful searches for pulsations from 
RX J1836. 2+5925 in the period range from 1 ms to 10 s and 
a pulsed-fraction upper limit of 35%. 

There are fewer observations of the GeV J1837-0610 
region. A 96-ms radio pulsar (PSR J1837-0604) was discov- 
ered and suggested to be associated wit h GeV J1837-0610 
(3EG J1837-0606) (|D'Amico et al.ll200lT ). However, the po- 
sition of PSR J1837-0604 is at the rim of the 95% confi- 
dence contour of GeV J1837-0610 and is inconsistent with 
AX J1837.5-0610, which is the single X-ray point source 
within the 95% confidence contour of GeV J1837-0610 in 
ASCA images. 

Chandra obs ervations of AX J2021. 1+3651 
l|Hessels et all |2004 ) indicate that this X-ray source is 
a pulsar wind nebula and there is an X-ray point source 
located at the center of this nebula, where a 104-ms radi o 
pulsar (PSR J2021+3651) was found l|Roberts et all 120021 ). 
The X-ray pulsation of this point source at the radio period 
was only marginally found in the Cha ndra data with an 
#-test score of 11.7 (|Hessels et al.ll2004l ). 

All the data were processed with standard procedures 
as outlined below. Photons for timing analysis were selected 
within a circle around the target and within a certain energy 
range. The circle radii and energy-range boundaries were so 
chosen that the target contribution was clearly detected. 
The target positions, the radii of the circles and the energy 
ranges adopted in this work are listed in Table [1] All the 
photon arrival times were then corrected to the solar system 
barycenter. The corrected time list was then used for timing 
analysis to look for periodicity. 
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Table 1. List of data employed in this study. Column 1 is the name of the X-ray source with its possible gamma-ray counterpart in 
the parentheses. Column 2 is the instrument used to obtain that data. Column 3 is the data epoch (in MJD), which is the middle point 
in the whole data time span. Column 4 and 5 are the total exposure time and the total time span of the data used. Both are in units 
of ks. Column 6 and 7 are the target position, photon numbers after extraction, the radius of the extraction circle, and the extraction 
energy range (in keV) for the X-ray target under study. Column 8 is the sequence number to label the data used in this work for later 
referencing. 
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3 THE CROSS-CHECKING METHOD FOR 
PERIODICITY SEARCH 

Periodicity search in sparse data usually employs photon- 
counting coherent-timing methods, that is, each photon ar- 
rival time is converted to a phase value with a given period 
and its time derivatives. The distribution of these phase val- 
ues is then tested against the null hypothesis with a certain 
statistical method to assess the probability that the observed 
distribution is consistent with a statisticall y flat distribu- 
tion. In this work, w e adopted the H-test ( De Jager et al.l 
Il989l ; iDe Jagerlll994h as the method to assess the random 
probability. 

The difficulty to find the period in sparse data, in ad- 
dition to the limited number of observations, is that a blind 
search with a huge number of trial periods (and period time 
derivatives) is unavoidable if an a priori period or period 
range is not known. When the total number of independent 
trials is taken into account, one seldom obtains a signifi- 
cantly low random probability. A possible way to improve 
upon this is to cross-check two different data to look for pos- 
sible support for the reality of a tempting candidate period 
found in one data set. 



For each of the X-ray sources discussed in the previous 
section we first performed a blind search in every available 
data of that source. Depending on the time resolution of 
the data different search ranges of trial periods were chosen. 
For all the ASCA/GIS data, the searching range is from 
0.1 s to 100 s. For all the others, it is from 0.01 s to 100 
s. Since each data set spans a relatively short time, we con- 
strained ourselves at this stage to consider trial periods only 
without their time derivatives. This makes the search a one- 
dimensional one. We chose the searching step to be a small 
fraction (typically 1/5) of the corresponding independent 
Fourier spacing, which is P 2 /T in the period domain (1/T 
in frequency), where P is the trial period and T is the total 
time span of the data. A blind search usually gives several 
tempting features. One example is shown in Fig.[T] For read- 
ers' convenience, the random probability for a single trial to 
have an H-vahie larger than a certain number is plotted in 
Fig.H 

From the result of a blind search in a certain data set, 
a few trial periods with high -ff-values were picked out for 
cross-checking. The number of such trial periods to pick out 
is obviously arbitrary. We picked out the top five in this 
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Figure 1. Examples of //-test blind search results. Shown here 
are those of AX J1420.1-6049 with data sets 11 and 12 (see Ta- 
ble [TJ. The number of independent trial frequencies is about 10 6 
for each case. Most of the trials yield a very low //-value. Only 
//-values larger than 20 are plotted. 
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Figure 2. The //-test probability distribution (De Jager et al.l 
Il989h . The plotted curve is the null hypothesis random probabil- 
ity for a single trial to have an //-value larger than Hq . 

work. Each of these five usually has a significantly low ran- 
dom probability if the total number of independent trials in 
that blind search is not taken into account. As mentioned 
above, a blind search has a huge number of independent tri- 
als, typically more than 10 6 . Therefore, for each of them, 
we look for possible associated features in other data sets 
to assess its reality. To do this, we further set two criteria. 
At first we look for associated features of the //-test results 
in another data set only within the range that gives a cor- 
responding characteristic age (jp) larger than 1000 years. 
A neutron star younger than that age probably would have 
a record in human history. In determining that range, the 
uncertainty of the p e riod i s estimated based on the method 
described in iLeahvl (1 19871 ). Secondly, the most prominent 
features in that range are identified and the resultant ran- 
dom probability is calculated, which is the product of the 
random probability of the feature in the first data set and 
that of the feature in the second data set. When calculating 
the random probability of the feature in the second data set, 



Figure 3. The identification of the candidate period at 159 ms 
for AX J1418. 6-6045 using data sets 04 and 05. The trial pe- 
riod range shown in the lower panel of this figure (from data set 
05) is roughly the range specified by the criterion that the corre- 
sponding characteristic age is larger than 1000 years if a possible 
period is to be associated with the one in the upper panel. The 
most significant feature in that range is at 0.15937378 sec, and 
its association with the one at 0.15935800 sec in the upper panel 
yields a resultant random probability of 5.3 X 10 — 3 . We list all 
such cases with a random probability lower than 0.01 as a candi- 
date period in Table [2] 

only independent trials within the range specified by the first 
criterion are considered. Oversampling was not taken into 
account in our computation of the random probability. We 
keep only those cases whose resultant random probability 
is less than 0.01. The periods and the corresponding period 
time derivatives for those cases are listed in Table [2] The 
probability threshold of 0.01 is arbitrary. Of course a higher 
one gives more candidate periods, but with less significance. 
Two examples of such an association are shown in Fig. [3] 
andg] 

In data set 07, the candidate period of AX J1418.7- 
6058 at 0.10171064 s has an //-value of 56.7, which has a 
very low random probability for a single trial. With our pro- 
cedure, there will be candidate periods satisfying our criteria 
in all the data sets 08, 09, and 10. Here we only report the 
most significant one. Its random probability listed in Ta- 
ble [5] was de rived using an //-valu e equal to 50 instead of 
56.7, since in lDe Jager et ail (|l989l ) the signle trial random 
probability was provided only up to H equal to 50. Two sim- 
ilar cases are the candidate period at 0.014799899 s for AX 
J1420. 1-6049, which has an //-value of 53.8 in data set 14, 
and that at 0.13247455 s for AX J1809. 8-2332, which has an 
//-value of 59.5 in data set 15. Although in our procedure 
only the most prominent feature in the cross-checking win- 
dow is picked out for association, for the candidate period 
of AX J1418.7-6058 at 0.10139550 sec in data set 07, the 
association with the period of the second highest //-value 
in the cross-checking window of data set 09 is also reported 
here, since it is very close to the candidate period in data 
set 07 and its //-value is not too different from the highest. 
That association is shown in Figure [5] The two candidate 
periods of AX J1809.8-2332 listed in Table [2] are related 
to each other. Based on their pulse profiles we suggest the 
period at 0.13247455 s be the second harmonic of the one 
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Table 2. Candidate periods. Column 1 is the candidate period (in seconds) obtained from cross-checking the two dat a sets l i sted i n 
column 7. The number in the parentheses is the l-cr uncertainty to the last digit, estimated with the method described in iLeahvl dl98Tt) . 
Those followed by a triangle (A) are candidate periods favoured by examing the consistency with properties of currently known gamma- 
ray pulsars, as discussed in Section 4. The epoch (MJD) of this period is listed in column 5, which is the epoch of the data set containing 
the stronger signature for pulsation in the two data sets. Column 2, 3 and 4 are the corresponding period time derivative (in units of 
10 — 13 s s — 1 ), characteristic age (in years) and spin-down power (in units of 10 36 erg s _1 ). Column 6 is the resultant random probability 
as explained in the text. 
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4.2x10- 


3 


07. 


10 


0.1443257(1) 


18.35(2) 


1250 


24.13(2) 


51223.68496 


1.9x10- 


■4 


07. 


08 


0.3274125(3) 


31.29(8) 


1650 


3.515(9) 


51223.68496 


9.5x10- 


3 


08. 


09 


AX J1420.1-6049 


0.014799899(6) 


0.392(3) 


6000 


478(4) 


52534.12055 


< 6.2 x 1Q- 3 


13, 


14 


0.019576788(3) 


2.616(3) 


1200 


1376(1) 


52289.79017 


7.2x10- 


3 


13. 


Id 


0.033549387(8) 


0.133(7) 


40000 


13.9(8) 


52289.79017 


1.2x10- 


3 


13. 


14 


0.14526767(4) 


4.84(2) 


4750 


6.24(2) 


51223.68495 


4.1x10- 


■5 


11. 


12 


0.14526767(4) A 


1.35(3) 


17000 


1.74(3) 


51223.68495 


5.0x10- 


3 


12. 


13 


0.1881326(4) 


10.1(3) 


2900 


6.0(2) 


52534.12055 


3.8x10- 


3 


13, 


11 


0.2540350(5) 


13.17(7) 


3100 


3.17(2) 


50317.00729 


5.0x10- 


3 


11. 


12 


AX J1809.8-2332 


0.13247455(3) 


16.439(5) 


1300 


27.792(8) 


50526.78845 


< 1.9 x 10~ 2 


15, 


17 


0.2649491(2) 


27.4(1) 


1500 


5.80(2) 


50526.78845 


7.4x10" 


4 


15. 


16 


AX J1836. 2+5928 


0.022425706(5) 


0.77(9) 


5000 


2.7(3) xlO 2 


53440.60818 


6.5x10- 


3 


19, 


21 


0.03350910(2) 


2(2) 


2000 


2(2) xlO 2 


53438.80937 


5.3x10- 


3 


20. 


21 


0.034677614(1) 


0.0201(4) 


270000 


1.90(4) 


50800.00499 


5.5x10- 


3 


18, 


19 


0.1679234(4) A 


0.16(2) 


170000 


0.13(1) 


53440.60818 


9.1x10- 


3 


18, 


21 


AX J1837.5-0610 


0.1736610(3) 


18.6(5) 


1500 


14.0(3) 


51105.10478 


2.8x10- 


3 


22, 


23 


AX J2021. 1+3651 



none 



at 0.264949 s if they are real. In such a case, the derived 
period derivative, characteristic age and spin-down power 
for the association between periodicity signatures found in 
data sets 15 and 17 near 0.1325 s will be 2, 1 and 0.25 times 
those listed in Table[5]respectively. It makes this association 
even more unfavoured if one compares its spin-down power 
with that of the Crab pulsar; see discussion in the next sec- 
tion. On the other hand, the strong periodicity signature 



(Figure [6J makes the candidate period at 0.264949 s very 
tempting. 

To examine whether the procedure described above is 
useful, we applied it to the case of RX J1856. 5-3754, whose 
X-ray emission is believed to be thermal emission from the 
surface of a neutron star. There were several XMM obser- 
vations of that source, but no periodicity was found until a 
period of 7.0 55 s with a random probabil ity of 6 x 10~ 4 was 
reported by iTiengo fc MereghettH l|2007l ). That result was 
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obtained by analysing the XMM data taken in October 2006. 
Applying our cross-checking method as described above to 
previous data, we found that the 7-s period was also picked 
out and suggested as a candidate period. It demonstrates, 
although not guarantees, the usefulness of our cross-checking 
method. 



50 E 
40 h 
30 f- 
20 h 
10 |,J 

o ili 



0.14526767(4) sec 



0.14530 

Period (sec) 



Figure 4. The identification of the candidate period at 145 ms 
for AX J1420. 1-6049 using data sets 11 and 12. The same as in 
Figure [3] the trial period range shown in this figure is roughly 
the range specified by the criterion that the characteristic age is 
larger than 1000 years. This figure is in fact a zoom-in view of 
Figure [l] around 6.88564 Hz but plotted in the period domain. 
This example, with a resultant random probability of 4.1 X 10~ 5 , 
is the most significant signature found in this study. 




Figure 5. The identification of the candidate period at 101.4 
ms for AX J1418. 7-6058 using data sets 07 and 09. The same 
as in Figure \3\ and [4] the trial period range shown in the lower 
panel of this figure is roughly the range specified by the criterion 
that the characteristic age is larger than 1000 years. The two 
possible period associations, which result in different period time 
derivatives, are both included in Table [2] 




Frequency (Hz) 



Figure 6. The ff-test blind search results of AX J1809.8-2332 
with data set 15. The two highest features correspond to periods 
at 0.13247455 s and 0.2649491 s respectively. 



4 DISCUSSION 

We applied the periodicity cross-checking procedure de- 
scribed in Section 3 to the 8 X-ray sources listed in Table Q] 
Candidate periods obtained with this method are listed in 
Table Future periodicity searches at X-rays or gamma- 
rays may be conducted around these candidate periods, in- 
stead of performing blind searches. Among all these candi- 
date periods, the one with the lowest random probability 
(4.1 xlO" 5 ) is 145.2297 ms for AX J1420. 1-6049. Its implied 
characteristic age, 4800 years, is between that of the Crab 
pulsar and of the Vela pulsar. However, its corresponding 
spin-down power is similar to that of the Vela pulsar, which 
in turn is about 100 times smaller than that of the Crab. 
If there is a correlation between the characteristic age and 
the spin-down power, this candidate period is not favoured 
because its associated spin-down power is too small for its 
4800-year characteristic age. 

Indeed, among the currently known gamma-ray pulsars, 
the spin-down power decreases with increasing characteris- 
tic ages. The Crab pulsar, with a characteristic age of 10 3 ' 1 
years, has a spin-down power about 10 38 ' 7 erg/s. That of the 
Vela pulsar are 10 4 ' 1 years and 10 36 S erg/s and of Geminga 
are 10 5 ' 5 years and 10 34 ' 5 erg/s. One possible way to fur- 
ther assess the reality of these candidate periods, therefore, 
is to compare their corresponding properties, such as the 
characteristic age and the spin-down power, with that of 
known gamma-ray pulsars. Along this line of thinking, the 
candidate periods at 203 ms for AX J1418.2-6047, at 159 
ms for AX J1418.6-6045, at 19.6 ms and 145 ms for AX 
J1420. 1-6049, and at 22.4 ms, 33.5 ms and 168 ms for AX 
J1836. 2+5928 are favoured, since their corresponding char- 
acteristic ages and spin-down powers are roughly consistent 
with that of known gamma-ray pulsars. 

Another property to consider i s the X-ray to gamma - 
ray energy spectral index, defined in lRoberte et all <|2001ah . 
which assumes an X-ray counterpart to those unidentified 
GeV sources and leads to three X-ray brightness categories 
for the GeV sources under study. In that scheme, the Crab 
pulsar is X-ray bright, the Vela pulsar is X-ray moderate, 
and Geminga is X-ray faint. The GeV sources associated 
with the aforementioned X-ray sources are all X-ray moder- 
ate (similar to Vela), except for GeV J1835+5921, which is 
X-ray faint (similar to Geminga). Examing the characteris- 
tic ages, we note that the four candidate periods at 203 ms 
for AX J1418. 2-6047, at 159 ms for AX J1418.6-6045, at 145 
ms for AX J1420.1-6049 and at 168 ms for AX J1836.2+5928 
survive the comparison. The candi date period at 127.5 ms 
for RX J0007.0+7302 reported by iLin fc Chanel <|2005h is 
also supported by these comparisons. 

These five X-ray sources are possible counterparts to 
three (or four) GeV sources, namely, GeV J0008+7304, 
GeV J1417-6100 (or a probable separate source 3EG J1420- 
6038), and GeV J1835+5921. For the former two (or three), 
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if the candidate periods suggested here turn out to be true, 
they will be more Vela-like, in terms of timing properties. It 
would suggest that the major cause for a gamma-ray pulsar 
to be radio quiet is likely geometrical, rather than intrinsic, 
and the gamma-ray beam is wider than the radio one so 
that for a certain viewing angle both the radio and gamma- 
ray emissions are observed, as for the case of Vela, and in 
another viewing angle only the gamma-ray emission is ob- 
served. Instead, if all the radio-quiet gamma-ray pulsars are 
similar to Geminga in their timing properties, the cause of 
being radio quiet, although not yet clear, is probably intrin- 
sic. 

Timing properties of not-yet-discovered gamma-ray pul- 
sars may not be similar to the known ones, particularly be- 
cause the number of the latter is still small. All the candidate 
periods listed in Table [2] may therefore still be helpful for 
pinning down the periodicity of emissions from these sources 
in the future when new X-ray data or gamma-ray data ob- 
tained with AGILE and GLAST become available. A larger 
sample of gamma-ray pulsars, either radio loud or quiet, 
will help us to understand better the radiation mechanisms 
in the magnetospheres of gamma-ray pulsars. 
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